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Enterovirus (EV) is one of the major causative agents of hand, foot, and mouth disease in the Pacific-Asia region. In particular, 
EV71 causes severe central nervous system infections, and the fatality rates from EV71 infection are high. Moreover, an outbreak 
of respiratory illnesses caused by an emerging EV, EV68, recently occurred among over 1,000 young children in the United States 
and was also associated with neurological infections. Although enterovirus has emerged as a considerable global public health 
threat, no antiviral drug for clinical use is available. In the present work, we screened our compound library for agents targeting 
viral protease and identified a peptidyl aldehyde, NK-1.8k, that inhibits the proliferation of different EV71 strains and one EV68 
strain and that had a 50% effective concentration of 90 nM. Low cytotoxicity (50% cytotoxic concentration, >200 pM) indicated 
a high selective index of over 2,000. We further characterized a single amino acid substitution inside protease 3C (3C?"°), N69S, 
which conferred EV71 resistance to NK-1.8k, possibly by increasing the flexibility of the substrate binding pocket of 3CP™°. The 
combination of NK-1.8k and an EV71 RNA-dependent RNA polymerase inhibitor or entry inhibitor exhibited a strong synergis- 
tic anti-EV71 effect. Our findings suggest that NK-1.8k could potentially be developed for anti-EV therapy. 


nterovirus (EV) and coxsackievirus (CV), which belong to the EV71 3CP"° (6, 20, 21), and pleconaril, GPP3, and cyclosporine 
Picornaviridae family, are the major causative agents of hand, block the entry of enteroviruses and coxsackieviruses (22-24). 
foot, and mouth disease (HFMD) in the Pacific-Asia region (1). | However, none of these reported inhibitors have been able to be 
Young children are more susceptible to enterovirus and coxsacki- advanced to clinical therapeutics. 
evirus infection (2, 3), but a recent case study reported that cox- In this work, we screened our library containing a number 
sackievirus also causes severe infectious diseases in adults (4). In of chemical compounds targeting viral proteases and identified 
particular, EV71 causes severe encephalitis, aseptic meningitis, 4 peptidyl aldehyde, NK-1.8k, that inhibited the proliferation 
pulmonary edema, acute flaccid paralysis, and myocarditis, which 6 f various strains of EV71, as well as one EV68 strain, in human 
lead to high fatality rates (5, 6). Moreover, in August 2014, an (RD and 293T) and monkey (Vero) cell lines and that had the 
emerging EV, EV68, caused mild to severe respiratory illnesses best 50% effective concentration (ECs) of 90 nM among the 
among 1,116 young children 1m the United States and was also compounds tested. We characterized the profile of EV71 resis- 
associated with neurological infections (7, 8). . tance to the compound and demonstrated that the amino acid 
Both enterovirus and coxsackievirus belong to the Enterovirus a j . i 
substitution conferring resistance occurs in the viral protease. 
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mately 7,400 nucleotides (10). Translation of the single open read- . i ; 
ing frame is initiated by ribosomes that use an internal ribosomal this compound could potentially be developed. for Ant EY 


entry site located in the 5’ untranslated region (UTR) region of the therapy. 

viral genome and gives rise to a polyprotein of approximately 250 

kDa (11-14). This polyprotein is further processed into four 

structural proteins (VP1 to VP4) to form the viral capsid and 

seven structural proteins (2A to 2C and 3A to 3D) by two viral 

proteases (2A and protease 3C [3CP"°] J together with a protease Received 9 January 2015 Returned for modification 9 February 2015 

precursor, 3CD (15). Accepted 13 February 2015 
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MATERIALS AND METHODS 


Cell lines and viruses. African green monkey kidney (Vero) cells (ATCC 
CCL-81), human rhabdomyosarcoma (RD) cells, and human embryonic 
kidney 293T cells were grown in Dulbecco’s modified Eagle’s medium 
(DMEM; Gibco) supplemented with 10% fetal bovine serum (FBS; 
Gibco) at 37°C in a humidified incubator with 5% CO}. 

Plasmids containing EV71 strain SK-EV006 (GenBank accession 
number AB469182.1) and EV71 strain SK-EV006 with a green fluorescent 
protein (GFP) reporter gene inserted into the genome (EV71-GFP) were 
donated by Satoshi Koike (Tokyo Metropolitan Institute of Medical Sci- 
ence) for the initial phenotype screening. A single-round reporter system 
[the EV71(FY)-Luc pseudotype virus system] with pseudotype EV71 con- 
taining plasmids pcDNA6-FY-capsid and pEV71-Luc-replicon lacking 
the P1 region was kindly supplied by Wenhui Li from the National Insti- 
tute of Biological Sciences, Beijing, China. Plasmids containing human 
EV71 strains AnHuil (GenBank accession number GQ994988.1) and 
BrCr (GenBank accession number U22521) were kindly provided by Bo 
Zhang and Hualin Wang, respectively, from the Wuhan Institute of Vi- 
rology. The complete genome of EV68 prototypic strain Fermon CA62-1 
(GenBank accession number AY426531) was chemically synthesized and 
was used to rescue living virus. 

Inhibitors and reagents. Chemical compounds in our screening li- 
brary and NITD008 (25) and GPP3 (23), which were used for the synergy 
studies, were synthesized. Stock solutions were prepared at a concentra- 
tion of 10 mM in distilled water with dimethyl sulfoxide (DMSO) (final 
concentration, 0.1% [vol/vol]) and before use were diluted to different 
concentrations with DMEM containing 10% FBS. 

A MEGAscript T7 high-yield transcription kit was purchased from 
Ambion. The TRIzol reagent and a SuperScript III first-strand synthesis 
system for reverse transcription-PCR (RT-PCR) were purchased from 
Invitrogen. A QuantiTect SYBR green RT-PCR kit was purchased from 
Qiagen. A cell viability and proliferation assay (WST-1) was purchased 
from Roche. 

Virus titration. For EV71-GFP, virus titers were determined by mea- 
suring the level of GFP expression in RD cells infected with EV71-GFP. 
Briefly, the measurement was performed by seeding 3 X 10* RD cells per 
well in 96-well microtiter plates. After overnight culture, EV71-GFP was 
serially diluted 10-fold with DMEM containing 10% FBS (107 '- to 1078- 
fold dilutions) and added to RD cells. The plates were then incubated at 
37°C in 5% CO). After 3 to 4 days, the level of GFP expression was mon- 
itored using an epifluorescence microscope. 

For living virus, virus titers were determined by endpoint dilution 
assays (EPDAs). Briefly, the titers were measured by seeding 3 X 10* RD 
cells per well in 96-well microtiter plates. After overnight culture, virus 
was serially diluted 10-fold with DMEM containing 10% FBS (10° '- to 
10 °-fold dilutions) and added to the RD cells. The plates were then 
incubated at 37°C in 5% CO}. Virus titers were determined after 3 to 4 
days. The virus titer, expressed as the 50% tissue culture infective dose 
(TCID;,), was determined using EPDAs. 

Preparation of EV71 pseudotype virus. The EV71 pseudotype virus 
[EV71(FY)-Luc] was produced by using the method developed as previ- 
ously described (26). Briefly, the plasmid pEV71-Luc subgenomic repli- 
con was linearized through digestion with the Sall restriction enzyme and 
was used as the template for RNA transcription. The EV71 replicon RNA 
transcripts were prepared in vitro using Promega MEGAscript kits. The 
pcDNA6-FY-capsid plasmid was transfected into 293T cells at 60 to 80% 
confluence. After 24 h posttransfection (hpi), EV71 subgenomic replicon 
RNA was then transfected using the Lipofectamine 2000 reagent (Invitro- 
gen). The EV71 pseudotype virus was harvested at 24 h after RNA trans- 
fection with two freeze-thaw cycles. To quantify the EV71 pseudotype 
virus, the stocks were serially diluted 10-fold and incubated with RD cells 
for 24 h at 37°C. The cells were then harvested, and the luminescence was 
detected with a Bright-Glo luciferase (Luc) assay system according to the 
manufacturer’s protocol (Promega). For the inhibition assay, the EV71 
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pseudotype virus was diluted to give a final concentration of 10° relative 
luminescence units (RLU) per well in a 96-well plate. 

Phenotype screening. RD cells were seeded at 3 X 10* cells per well in 
96-well microtiter plates. After overnight culture at 37°C in 5% CO,, RD 
cells were treated with a gradient of concentrations of different com- 
pounds ranging from 0.045 to 100 pM. After 2 h, RD cells were infected 
with EV71-GFP at a multiplicity of infection (MOI) of 1. At 24 hpi, the 
level of GFP expression was monitored using an epifluorescence micro- 
scope. 

Western blot analysis. RD cells were treated with NK-1.8k at concen- 
trations ranging from 0.05 to 0.2 pM and lysed in a buffer containing (50 
mM Tris-HCl, pH 8.0, 150 mM NaCl, 1% NP-40, 0.5% sodium deoxy- 
cholate, 0.1% sodium dodecyl sulfate [SDS], 2 mM EDTA, 1 mM NaVO,, 
10 mM NaF). The level of expression of the EV71 VP1 protein in the 
lysates was determined by a spectrophotometer. Proteins were separated 
by SDS-polyacrylamide gel electrophoresis (PAGE) and transferred to a 
nitrocellulose membrane (Millipore). These membranes were blocked 
with a 5% nonfat dry milk solution in Tris-buffered saline for 4 h, incu- 
bated with specific primary antibodies (catalog number ab36367; Ab- 
cam), and then subsequently incubated with secondary antibodies conju- 
gated with horseradish peroxidase (HRP). Proteins were visualized by 
chemiluminescence using Clarity Western ECL substrate (Bio-Rad). 

qRT-PCR. The reduction in the amount of the viral genome in host 
cells after treatment with one of the test compounds was determined using 
living EV71 and cells in a quantitative RT-PCR (qRT-PCR)-based assay. 
Briefly, 10° cells were seeded in each well of 24-well tissue culture plates 
with complete culture medium, and the plates were incubated at 37°C in 
5% CO,. The medium was then replaced with medium containing EV71 
strain A, B, or C, as well as EV68 strain Fermon CA62-1, at an MOI of 1 
and compounds serially diluted to concentrations ranging from 0.0039 to 
2 pM in the presence of 10% FBS and 0.5% DMSO, respectively. After the 
cells were treated for 24 h, total cellular RNA was extracted with the TRI- 
zol reagent using standard protocols. A quantitative RT-PCR assay (for 
GADPH [glyceraldehyde-3-phosphate dehydrogenase], the primer se- 
quences were 5'-CCCACTCCTCCACCTTTGACG-3’ [forward primer] 
and 5'-CACCACCCTGTTGCTGTAGCCA-3’ [reverse primer] ) was per- 
formed using a QuantiTect SYBR green RT-PCR kit (Qiagen, Valencia, 
CA) following the manufacturer’s protocol. EVs and GADPH transcript 
levels were determined by the AAC; threshold cycle (C+) method. The 
percent inhibition was calculated as follows: [1 — (average number of 
compound-treated cells)/(average number of control cells)] X 100. Each 
data point represents the average of three replicates in cell culture. 

Inhibitory effect on different cells. In order to eliminate the differ- 
ences between cells, the EC;9s for RD, Vero, and 293T cells were deter- 
mined. RD (3 X 104 per well), Vero (2 X 104 per well), and 293T (3 X 104 
per well) cells were seeded in 96-well plates and cultured at 37°C in 5% 
CO, overnight. Cells were treated with a gradient of dilutions of the dif- 
ferent compounds ranging in concentration from 0.003 to 100 uM. EV71 
Fuyang was added to 96-well plates at an MOI of 1. The RNA was detected 
by qRT-PCR after 24 hpi, and the EC.,,s on different cells were calculated 
by the use of GraphPad Prism software. 

Virus titers in supernatants. The levels of reduction of the virus titers 
in the supernatants were determined by EPDAs. Briefly, RD cells (3 X 10° 
per well) were cultured in a 24-well plate and incubated at 37°C under 5% 
CO,. On the following days, RD cells were treated with serially diluted 
compounds at concentrations ranging from 0.0039 to 2 uM in the pres- 
ence of 10% FBS and 0.5% DMSO. After 2 h, EV71 (Fuyang strain) at an 
MOI of 1 was used to infect the cells. The viruses in the supernatants of the 
wells containing each concentration were added to new 96-well plates and 
diluted from 10! to 1078. The cytopathic effect was observed in every 
well after 72 hpi. The TCID.,s were calculated by use of the Reed and 
Muench calculator. 

Cytotoxicity. Cytotoxicity was measured by a WST-1-based cell pro- 
liferation assay using a cytotoxicity assay kit according to the manufactur- 
er’s protocol (catalog no. 05015944001; Roche, USA). Briefly, Vero cells 
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FIG 1 NK-1.8k inhibited EV71 proliferation. (A) Chemical structures of NK-1.8k and rupintrivir. (B) Concentration-dependent reduction of EV71 prolifera- 
tion following treatment with NK-1.8k. RD cells were infected with EV71-GFP at an MOI of 1 and treated with various concentrations of NK-1.8k (0.045 to 0.39 
uM) for 24 h or left untreated. DAPI (4’,6-diamidino-2-phenylindole) was used to visualize the nuclei (top) of the host cells, and enhanced GFP was used to 
monitor virus growth (bottom). (C) Inhibitory effect of NK-1.8k on EV71 RNA synthesis measured by qRT-PCR. RD cells were infected with EV71 at an MOI 
of 1 and treated with NK-1.8k at various concentrations (0.0039 to 2 uM) for 24 h or left untreated. Each point presents the average results from three 
independent experiments, and error bars represent SEMs (n = 3). (D) Dose-dependent reduction of EV71 VP1 expression. The concentration of NK-1.8k in the 
four lanes was 0.05, 0.1, 0.15, and 0.2 uM, from left to right, respectively. The level of GADPH expression was not affected by treatment with NK-1.8k. 


(2 X 10* per well) and 293T and RD cells (3 X 10* per well in 100 pl 10% 
FBS-DMEM) were cultured at 37°C under 5% CO, in 96-well plates, 
followed by the addition of 50 wl of a solution containing each of the 
compounds to be tested at a concentration ranging from 0.78 to 200 uM. 
The cells were incubated at 37°C for 24 h, and then 10 ul WST-1 reagent 
was added to each well. The luminescent signals were read at 490 nm with 
a microplate reader (Bio-Rad, USA). The percentage of viable cells after 
treatment with different concentrations of compound was calculated as 
follows: percentage of viable cells = 100 X (OD joo of the treated sample/ 
OD joo of the cell control sample), where OD 499 is the optical density 
(OD) at 490 nm. 

Time-of-addition assay. The time-of-addition effect was examined 
for NK-1.8k, NITD008, and GPP3. RD cells (3 X 10* cells per well in 100 
ul 10% FBS-DMEM) were cultured overnight at 37°C under 5% CO, in 
96-well plates. Subsequently, the cells were treated with 0.2 uM NK-1.8k, 
2 uM NITD008, and 1 pM GPP3 either concurrently with 100 TCID.9s of 
EV71(FY)-Luc (0 h) or at intervals of —6, —4, —2, 0, 2, 4, 6, 8, and 10 hpi. 
Another 96-well plate was treated with NK-1.8k (0.2 uM), NITD008 (2 
uM), and GPP3 (1 pM) and infected with EV71(FY)-Luc after 6 h. After 
incubation at 37°C for 24 h, antiviral activity was determined by measure- 
ment of the reduction of luciferase activity compared with the luciferase 
activity of the control cultures. 

Selection of drug-resistant virus. RD cells were seeded at 6 X 10° 
cells/well in 6-well plates. On the following day, the medium was removed 
and replaced with DMEM containing 10% FBS and 0.2 uM NK-1.8k; 
0.5% DMSO was used as a control. After 4 h, EV71 (Fuyang strain) was 
used to infect RD cells at an MOI of 0.1. After 3 days, the virus (termed 
passage 1 virus) was harvested when RD cells reached 70 to 90% conflu- 
ence. A fresh 96-well plate was prepared and replaced with 0.2 uM NK- 
1.8k. RD cells were infected with EV71 at an MOI of 0.1, and virus from 
passage 2 was harvested. After three rounds of passage, the concentration 
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of NK-1.8k was increased to 0.4 uM for another three rounds of passage 
and then to 0.8 uM for two rounds. 

Identification of drug resistance-conferring mutations. For analysis 
of EV71 RNA resistance-conferring mutations, cellular RNA was ex- 
tracted with the TRIzol reagent (Invitrogen) according to the manufac- 
turer’s instructions. For reverse transcription-PCR, first-strand cDNA 
was synthesized using a gene-specific primer (5’-ACCCCCACCAGTCA 
CATTCACG-3’) and the SuperScript III first-strand synthesis system for 
RT-PCR (Invitrogen) according to the manufacturer’s instructions. The 
protein-coding region of the EV71 genome was amplified by PCR in 5 
short fragments: fragment 1 was amplified with primers EV71-718-sense 
(5'-ATCTTGACCCTTAACACAGC-3') and EV71-2046-anti (5’-GACC 
ATTGGGTGTAGTACCC-3’), fragment 2 was amplified with primers 
EV71-1975-sense (5'-CGATCCTGGGCGAAGTGGAC-3’) and EV71- 
3345-anti (5’-TGITGTCCAAATT TCCCAAG-3’), fragment 3 was 
amplified with primers EV71-3248-sense (5'-TACCTATTCAAAGCCAA 
CCC-3') and EV71-4643-anti (5'-ATAAAGACATATCCTTGCCG-3’), 
fragment 4 was amplified with primers EV71-4569-sense (5’-ACGGCTA 
CAAGCAACAGGTG-3’) and EV71-6044-anti (5’-TTCCCTCGAAGAT 
ATCATGG-3’), and fragment 5 was amplified with primers EV71-5972- 
sense (5'-GGAAGGCTCAACATCAATGG-3’) and EV71-7345-anti (5'- 
GGGTTGAGGTGTGTATAGCC-3’). The RT-PCR products of the 5 
fragments of resistant EV71 or control EV71 were ligated into the TA 
cloning vector PMD 18-T (TaKaRa). For each fragment, multiple individ- 
ual bacterial colonies were isolated, and the purified plasmid DNA was 
sequenced. The sequences were aligned using Sequencher (version 5.0) 
and BioEdit software. 

Sensitivity of EV71 WT and N69S mutant virus to NK-1.8k. The 
recombinant plasmid carrying EV71 (strain Fuyang) with a T7 promoter, 
a HindIII site before the 5’ UTR, and an Xbal site in the poly(A) sequence 
after the 3’ end was ligated into pcDNA 3.1. The N69S amino acid substi- 
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FIG2 Anti-EV71 activities of NK-1.8k at different times of addition. (A to C) Inhibition of EV71(FY)-Luc pseudotype virus infection of RD cells by NK-1.8k (0.2 
uM) (A), NITD008 (2 uM) (B), and GPP3 (1 pM) (C) at various times of addition. The luciferase levels at 24 hpi were quantified by measuring firefly luciferase 
activity in relative luminescence units (RLU). (D) EV71 subgenomic replicon RNA lacking the P1 region was transfected into RD cells treated with NK-1.8k (0.2 
uM), NITDO008 (2 uM), and GPP3 (1 pM) or left untreated (control). The luciferase levels were quantified at 24 hpi. Each point represents the average results 


from three independent experiments, and error bars represent SEMs (n = 3). 


tution was generated by using a Fast mutagenesis system kit (Transgen) 
with forward primer 5’-ATGAGCAAGGAGTCAGCTTGGAATTAACC 
C-3’ and reverse primer 5’-CTGACTCCTTGCTCATCCACTAGCTCAA 
CTG-3’. The full-length RNAs of wild-type (WT) EV71 and the N69S 
mutant were transcribed from recombinant plasmid linearized with Xbal 
in vitro using a RiboMAX large-scale RNA production system (T7 kit; 
Promega). RD cells were electroporated with 10 wg of WT or mutant 
genome-length RNA. The WT and mutant viruses were harvested at 24 h 
posttransfection. RD cells (3 X 104 cells per well) were seeded in a 96-well 
plate and cultured at 37°C for 24 h. The WT and mutant virus titers were 
determined by EPDAs as described above. For testing the inhibitory ac- 
tivity of NK-1.8k between the WT and compound-resistant virus, 500 pl 
of the RD cell suspension was immediately seeded in a 24-well plate and 
the cells were incubated for 24 h. Cells were treated with serial dilutions of 
NK-1.8k and infected with WT or N69S virus. At 24 hpi, RNAs were 
extracted and quantify by RT-PCR. 

Interaction of NK-1.8k with other inhibitors of EV71 proliferation. 
The interaction of NK-1.8k and the other inhibitors was studied by use of 
a previously reported protocol and mathematical model, with some mod- 
ifications (27, 28). Briefly, RD cells were seeded at 3 X 104 cells per well in 
96-well microtiter plates. After an overnight culture at 37°C, a different 
final concentration of NK-1.8k (0, 0.3125, 0.0625, 0.125, 0.25, and 0.5 
uM) was added to each row of the 96-well plate. Simultaneously, a differ- 
ent final concentration of NITD008 (0, 0.0625, 0.125, 0.25, 0.5, 1, 2, and 4 
uM) was added to each column of the plate. RD cells were infected with 
the EV71 pseudotype virus after treatment with the two inhibitors and 
incubated for 24 h. The combination of NK-1.8k and GPP3 was tested by 
the same method, but with some adjustment. The concentrations of GPP3 
used were 0.0005, 0.0019, 0.0078, 0.0312, 0.125, 0.5, and 2 uM. Antiviral 
activities were determined by measuring the reduction of luciferase activ- 
ity in the cells. 

Differential surface plots at the 95% confidence interval (CI) were 
calculated and generated by using the MacSynergy II program for the 
drug-drug interaction according to the Bliss independence model (28). 
The combination’s effect was determined by subtracting the experimental 
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values from the theoretical additive values (27, 28). In the three-dimen- 
sional differential surface plot, the peaks above and below a theoretical 
plane for additivity represent synergy and antagonism, respectively (27, 
28). The data sets were interpreted as follows: volumes of synergy or an- 
tagonism with values of <25 mM?% are insignificant, those with values of 
25 to 50 mM?% are minor but significant, those with values of 50 to 100 
mM?% are moderate and probably important in vivo, and those with 
values of >100 mM*% are strong and likely to be important in vivo. 

Protein production. Measurement of the production of EV71 WT, 
N69S mutant, and C147S mutant 3CP"° was based on a previously re- 
ported protocol (29-32). Briefly, the genomic region that codes for EV71 
3C was amplified and cloned into the pET-28a vector (GE Healthcare) 
using the NdeI and NotI restriction sites. The recombinant plasmids were 
transformed into Escherichia coli strain BL21(DE3). Kanamycin-resistant 
colonies were grown in Luria-Bertani medium at 37°C until the OD¢q9 
reached 0.8. Isopropyl-B-p-1-thiogalactopyranoside was added to a final 
concentration of 0.5 mM, and the cultures were grown for an additional 
20 h at 18°C. Cells were harvested by centrifugation, resuspended, and 
homogenized in lysis buffer containing 50 mM Tris-HCl (pH 7.0) and 200 
mM NaCl, using a low-temperature ultra-high-pressure cell disrupter 
(JNBIO, Guangzhou, China). The lysate was centrifuged at 20,000 X g for 
30 min to remove the cell debris. The supernatant was loaded onto an 
Ni-nitrilotriacetic acid affinity chromatography column (Qiagen). After 
washing with 50 mM Tris-HCl (pH 7.0), 200 mM NaCl, and 50 mM 
imidazole, the 3CP™° protein was washed with 50 mM Tris-HCl (pH 7.0), 
200 mM NaCl, and 500 mM imidazole. The sample was further purified 
using a HiTrap S column (GE Healthcare) with a linear gradient of from 0 
mM to 500 mM NaCl with 50 mM Tris-HCl (pH 7.0) and 2 mM dithio- 
threitol (DTT). The target proteins were concentrated to 12 mg/ml in a 
buffer with 20 mM Tris-HCl (pH 7.0), 500 mM NaCl, and 2 mM DTT for 
storage. 

In vitro protease activity assay. In the in vitro cleavage assay, the 
peptide was attached to a fluorescence-quenching pair of peptides, and 
the fluorogenic peptide NMA-IEALFQGPPK(DNP)FR (where NMA is 
N-methylanthranyl and DNP is dinitrophenol) was used as the substrate 
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FIG 3 Quantification of anti-EV71 effect of NK-1.8k on different cell lines. (A, D, and G) RD, 293T, and Vero cells were treated with NK-1.8k and infected with 
EV71 strain Fuyang. EV71 RNA levels were measured by qRT-PCR at 24 hpi. (B, E, and H) RD, 293T, and Vero cells were treated with NK-1.8k and infected with 
EV71 strain Fuyang. The virus titers in the supernatants were determined at 72 hpi. (C, F, and I) Cytotoxicity in different cell lines. RD, 293T, and Vero cells were 
incubated with NK-1.8k at the indicated concentrations for 24 h. Cell viability was measured using a WST- 1-based assay. Each point represents the average results 


from three independent experiments, and error bars represent SEMs (n = 3). 


for the enzyme assay (33). The assays were performed with 100-1 samples 
containing 0.5 uM EV71 WT protease or the N69S or C147S mutant 
protease and 20 uM substrate in 50 mM Tris, pH 7.0, buffer. The enzymes 
were preincubated at 30°C for 30 min prior to the addition of the sub- 
strate. The reaction mixture was then incubated for 2 h at 30°C. The 
reaction was terminated by the addition of acetonitrile at 10-min inter- 
vals. The intensity of the OD for NMA and DNP was read at an excitation 
wavelength (A,,.) of 340 nm and an emission wavelength (Nom) of 440 nm. 


RESULTS 

NK-1.8k is an inhibitor of EV71 infection. Our aim in this work 
was to find an inhibitor that can inhibit the correct function of 
viral 3CP’° during the replication of EV71. A quick phenotype 
screening was performed using a total of 87 protease inhibitor 
derivatives that were initially assayed for their antiviral activity 
against EV71 by using RD cells infected with EV71 strain SK- 
EV006 that expressed GFP upon infection (34). Most of the tested 
compounds exhibited no obvious anti-EV71 effect at a concentra- 
tion of 20 uM, and GFP expression remained obvious in infected 
RD cells. In contrast, treatment with a peptidyl aldehyde, NK-1.8k 
{N-[(S)-3-(4-fluorophenyl)-1-oxo-1-{(S)-1-oxo0-3-[(S)-2- 
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oxopiperidin-3-yl] propan-2-ylamino}propan-2-yl]-5-methyl- 
isoxazole-3-carboxamide} (Fig. 1A), achieved apparent reduc- 
tions in the level of expression of the reporter GFP at a concentra- 
tion of 0.09 uM (Fig. 1B). NK-1.8k was further shown to cause 
concentration-dependent protection against EV71 infection by 
measurement of the level of GFP expression (Fig. 1B). 

NK-1.8k is a novel peptidomimetic compound which has 
structural characteristics that distinguish it from rupintrivir (21) 
(Fig. 1A). NK-1.8k features a six-member-ring lactam coupled 
with an aldehyde functional group, while rupintrivir has a five- 
member-ring lactam coupled with an a,f-unsaturated ester. 
Moreover, rupintrivir is considered a tripeptidic analogue, 
whereas NK-1.8k is similar to a dipeptide. These features indicate 
that NK-1.8k has better drug-like properties with a lower molec- 
ular weight and fewer hydrogen bond donors and acceptors than 
rupintrivir (35). 

We thus further studied the anti-EV71 activity of NK-1.8k by 
infecting RD cells with EV71 strain SK-EV006. The infected cells 
were treated with various concentrations of NK-1.8k for 24 h, and 
the remaining levels of EV71 RNA in the cells were determined by 
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FIG 4 Quantification of antiviral effect of NK-1.8k on different virus strains. RD cells were treated with NK-1.8k and infected with living EV71 strain BrCr 
(genotype A) (A), EV71 strain SK-EV006 (genotype B) (B), EV71 strain Fuyang (genotype C) (C), and EV68 strain Fermon CA62-1 (D). The remaining levels of 
viral RNA in the cells were determined by the qRT-PCR method. Each point represents the average results from three independent experiments, and error bars 


represent SEMs (n = 3). 


qRT-PCR. The results revealed that NK-1.8k showed a concentra- 
tion-dependent reduction of EV71 RNA in infected cells with an 
EC;, of 34.5 nM (Fig. 1C). Moreover, it was observed that the 
expression of the VP1 protein of EV71 (the major component of 
the viral capsid) but not that of host GADPH, a host housekeeping 
gene, in the infected cells was also decreased by treatment with 
NK-1.8k (Fig. 1D). All these results support the possibility that 
NK-1.8k has an inhibitory effect on EV71 proliferation. 

NK-1.8kimpacts EV71 genome replication. Because NK-1.8k 
is potentially a protease inhibitor, we first analyzed whether NK- 
1.8k affects viral replication but not the entry process. A previ- 
ously reported pseudotype EV71 system was used to exclude the 
possibility of an impact of virus reinfection (36). 

We infected RD cells with EV71(FY)-Luc pseudotype virus and 
then treated the cells with NK-1.8k at a concentration of 0.2 uM, 
NITD008 (a reported polymerase inhibitor) (18, 19) at a concen- 
tration of 4 mM, and GPP3 (23) at a concentration of 1 pM at —6, 
—4, —2, 2, 4, 6, 8, and 10 hpi, in which 0 hpi indicates the time that 
the inhibitors were supplied immediately after virus infection. The 
results showed that the inhibition of EV71 by NK-1.8k had no 
dependence on treatment time before 10 hpi (Fig. 2A). The treat- 
ments with NK-1.8k at —6 to 8 hpi displayed the strongest antivi- 
ral effect, which is similar to the findings reported for the viral 
replication inhibitor NITD008 (Fig. 2B). In sharp contrast, the 
antiviral effect of GPP3 showed a dramatic decrease from 2 hpi 
(Fig. 20). 

We next directly transfected EV71 subgenomic replicon RNA 
lacking the P1 region into RD cells treated with NK-1.8k (0.2 uM), 
NITD008 (2 uM), and GPP3 (1 pM). The results showed that 
NK-1.8k caused the inhibition of EV71 replication to a level sim- 
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ilar to the level of inhibition caused by NITD008 but distinctly 
different from that caused by GPP3 (Fig. 2D). Because this trans- 
fection excluded the region of the genome responsible for the 
entry step of viral infection, all of these results support the sugges- 
tion that NK-1.8k inhibits EV71 proliferation by affecting the viral 
replication stage. 

Quantification of the anti-EV71 effect of NK-1.8k in different 
cell lines. To assess possible cell type- and species-dependent dif- 
ferences in the anti-EV71 activity of NK-1.8k, its antiviral effects 
in RD, 293T, and Vero cells were examined (Fig. 3). 

The results revealed that NK-1.8k suppresses the viral RNA 
yield of EV71-infected RD, 293T, and Vero cells in a dose-depen- 
dent fashion. The EC.9s of NK-1.8k against EV71 SK-EV006 in- 
fection in RD, 293T, and Vero cells were determined to be 0.108, 
0.108 and 2.41 uM, respectively (Fig. 3A, D, and G). 

We also measured the reduction of virus titers after treatment 
with NK-1.8k in EV71-infected RD cells. Consistently, the results 
indicated that the virus titers in the supernatants were also signif- 
icantly attenuated by treatment with NK-1.8k (Fig. 3B, E, and H). 

The cytotoxicity of NK-1.8k for RD, 293T, and Vero cells was 
also tested. The results revealed that NK-1.8k at a high concentra- 
tion of 200 uM has no obvious toxicity for any of the three cell 
lines (Fig. 3C, F, and I). Taken together, these results indicate that 
NK-1.8k inhibited EV71 with a high potency and the anti-EV71 
activity of NK-1.8k has a high to satisfactory selectivity index (SI; 
where SI is equal to CC59/EC59, where CCo is the 50% cytotoxic 
concentration) of up to 1,000 in RD cells. 

NK-1.8k has a broad antiviral effect against multiple EV71 
strains and EV68. EV71 is generally characterized as three geno- 
types, i.e., genotypes A, B, and C. We hypothesized that NK-1.8k 
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FIG 5 Mutation identified in NK-1.8k-resistant EV71. (A) Mutations recov- 
ered from resistant viruses were tabulated according to the amino acid posi- 
tions and the corresponding protein. Sequencing of the complete genomes of 
24 viruses from SelA, SelB, and SelC was performed. (B) Sequence chromato- 
grams of the mutations found in resistant viruses. Representative sequencing 
chromatograms are presented for SelA, SelB, and SelC. 


may have a broad antiviral effect against viruses of all three geno- 
types. Therefore, we infected RD cells with three different virus 
strains, including BrCr (genotype A), SK-EV006 (genotype B), 
and Fuyang (genotype C), treated them with NK-1.8k, and mea- 
sured the remaining levels of EV71 RNA inside host cells through 
a qRT-PCR method. The results showed that NITD008 inhibited 
the viruses of strains BrCr, SK-EV006, and Fuyang with EC.s of 
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0.105, 0.093, and 0.102 uM, respectively (Fig. 4A to C). Moreover, 
we also tested the antiviral effect of NK-1.8k on EV68 strain Fer- 
mon CA62-1 and found an EC; of 0.151 uM against EV68 infec- 
tion. All these results suggest that NK-1.8k has a broad range of 
inhibition and can inhibit different enteroviruses. 

Selection and identification of resistant viruses. The NK- 
1.8k-resistant viruses were selected by continuously culturing 
EV71 for eight rounds in the presence of gradually increasing con- 
centration of NK-1.8k: 0.2 uM (three rounds), 0.4 uM (three 
rounds), and 0.8 uM (two rounds). Three independent selections 
(SelA, SelB, and SelC) were performed to test the reproducibility 
of the resistance results. 

Sequence analyses of the entire genomes of multiple resistant 
viruses identified a single A-to-G mutation at nucleotide position 
5592 of the whole EV71 genome. This mutation translated into a 
single amino acid replacement of an asparagine residue by a serine 
residue at position 69 within the 3CP"° polypeptide (Fig. 5A and 
B). Because NK-1.8k was screened out of the library of com- 
pounds targeting the viral protease, this drug resistance is reason- 
able. 

Sensitivity of WT EV71 and EV71 with the N69S mutation to 
NK-1.8k. We further examined the fitness and drug resistance of 
the N69S mutant virus. We first examined the effect of this muta- 
tion on viral replication in the absence of the compound. RD cells 
were transfected with equal amounts of WT EV71 or EV71 mutant 
N69S RNA and assayed for virus titers after transfection. The re- 
sults revealed that the titer of the N69S mutant virus had an ap- 
parent decrease in comparison to that of the WT virus, suggesting 
that the fitness of the N69S mutant virus is inferior to that of the 
WT virus (Fig. 6A). 

We next examined the sensitivity of mutant viruses to NK-1.8k 
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FIG 6 Sensitivity of WT EV71 and the N69S mutant virus to NK-1.8k. (A) Titers of EV71 WT and N69S mutant viruses. (B) RD cells were treated with NK-1.8k 
at concentrations ranging from 2 to 0.0039 uM and infected with the EV71 WT or N69S mutant at an MOI of 1. The inhibitory activities were determined by 
qRT-PCR. (C) Protease activity of EV71 WT, the N69S mutant, and the C147S mutant. The intensities of the ODs for EV71 WT, N69S mutant, and C147S mutant 
proteases were read at a Nex of 340 nm and a Nem of 440 nm every 10 min for 120 min. AU, absorbance units. 
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FIG7 Interaction of NK-1.8k with other inhibitors of EV71 proliferation. (A and B) Anti-EV71 effects of the combinations NK-1.8k plus NITD008 and NK-1.8k 
plus GPP3. (C and D) Differential surface plots at the 95% confidence interval (CI) were calculated and generated by using the MacSynergy II program for the 
interactions of the combinations NK-1.8k plus NITD008 (C) and NK-1.8k plus GPP3 (D). The y axes show the volumes of synergy or antagonism (in mM?%). 


inhibition using a qRT-PCR assay. The results showed that the 
N69S mutation conferred complete resistance to NK-1.8k, but the 
WT virus did not show any resistance (Fig. 6B). Collectively, all 
these results demonstrate that the N69S mutation attenuates the 
fitness of virus growth but helps the virus to escape inhibition by 
NK-1.8k. 

Protease activities of WT and N69S mutant 3C?". In order 
to verify the influence of the EV71 N69S amino acid substitution 
on enzyme activity, the peptide NMA-IEALFQGPPK(DNP)FR 
was used as the substrate for the protease assay. The OD intensity 
was read at an Nex of 340 nm and an Nem of 440 nm every 10 min. 

The results determined from the change in the intensity of the 
OD within 120 min of exposure demonstrated that the 3CP”® pro- 
tease with the N69S mutation retained a low level of protease 
activity compared to that of the catalytically defective mutant with 
the C147S mutation but that its protease activity was much lower 
than that of WT 3CP™ (Fig. 6C). This is consistent with the lower 
fitness of the virus with the N69S mutation. 

Effects of the combination of NK-1.8k and polymerase and 
entry inhibitors against EV71 infection. The use of therapies 
with combinations of drugs with different mechanisms of action 
has been proven to be an effective approach to the treatment of 
viral infections, including HIV and hepatitis C virus infections. 
We reasoned that the use of a combination of entry and polymer- 
ase inhibitors of EV71 with NK-1.8k might offer benefits over the 
use of each single monotherapy. To test this hypothesis, we mea- 
sured the inhibition of EV71 propagation by NITD008 and GPP3 
in combination with NK-1.8k. 

The combinations of NK-1.8k with NITD008 and with GPP3 
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were evaluated with the single-round pseudotype EV71 reporter 
system. The results revealed a concentration-dependent inhibi- 
tion of EV71 replication by NK-1.8k, NITD008, or GPP3 alone or 
with any two of these drugs in combination. The combinations of 
NK-1.8k plus NITD008 (Fig. 7A) and NK-1.8k plus GPP3 (Fig. 
7B) at various concentrations resulted in inhibition greater than 
that achieved with either compound in the combination alone. 

The data were further analyzed using a mathematical model, 
MacSynergy II, to determine whether the effect of the combina- 
tions was synergistic, additive, or antagonistic. The results re- 
vealed that the combination of NK-1.8k and NITD008 had anti- 
viral effects that were significantly more potent than the 
theoretical additive effects above concentrations of 0.25 uM for 
NITD008 and 1.9 pM for GPP3. This finding supports the sugges- 
tion that this combination is indeed synergistic (Fig. 7C and D). 
Because the effects of NK-1.8k, NITD008, and GPP3 on cell via- 
bility were previously evaluated to ensure that the inhibition of 
EV71 was not due to cytotoxicity, it is not surprising that no sig- 
nificant cytotoxicity was observed with the combination of NK- 
1.8k and either of the other two compounds (data not shown). 
Therefore, the observed synergy is indeed specific and does not 
reflect synergistic toxicity. 


DISCUSSION 


All picornaviruses, such as EV71, polioviruses, rhinoviruses 
(RVs), CVs, and hepatitis A virus, encode a cysteine protease 
(named 3C?"°) that proteolytically processes the polyprotein into 
functional viral proteins (37). Because the viral protease has strict 
substrate specificity and lacks homologues in mammalian cells, it 
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FIG 8 N69S may alter the conformation of the active site of EV71 3C?"°. The structures of WT EV71 3CP"® (PDB accession number 3OSY) (A), the catalytically 
defective C147A mutant of CVA16 3CP™° in complex with a substrate peptide (PDB accession number 3SJ9) (B), EV71 3C?*° in complex with rupintrivir 
(AG7088; PDB accession number 3SJO) (C), and a modeled structure of EV71 3CP"° with the N69S mutant (D) were aligned and are presented in the same 
orientation. The polypeptides of protease are shown as white cartoon diagrams; the substrate peptide and rupintrivir (AG7088) molecule are displayed as stick 
models. Key residues N69 (or N69S), E71, T132, and C147 are also represented as colored sticks. The hydrogen bonds between the side chains of N69 and 
E71/T132 are denoted as dashed lines. (E) Sequence alignment of EV71 strain Fuyang and EV68 strain Fermon CA62-1. The strictly conserved residues are shown 
as white text with a red background, whereas the conserved residues are displayed as red text with a white background. The blue frame highlights drug 


resistance-conferring residue N69. 


is recognized to be the most essential target, together with viral 
polymerase, for antiviral discovery (38). Therefore, we systemati- 
cally analyzed the substrate specificity of virus-encoded cysteine 
proteases and generated a chemical compound library targeting 
viral cysteine protease activity through a combinational chemistry 
method. 

In this study, by screening the compounds in our protease- 
targeted library, we verified that a peptidyl aldehyde, NK-1.8k, 
targets the EV71 protease to inhibit the replication of various virus 
strains. NK-1.8k has no obvious cytotoxicity, even at a high con- 
centration of 200 pM, for different cell lines. Through the screen- 
ing of drug-resistant viruses, we identified the amino acid substi- 
tution N69S on the EV71 protease to be the substitution that 
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confers resistance. Although this substitution attenuated the pro- 
tease activity of 3CP™ and the fitness of EV71, it helped the virus to 
completely escape from NK-1.8k treatment, supporting the sug- 
gestion that NK-1.8k directly inhibits the activity of viral protease 
itself. 

A previous study identified several mutations in human rhino- 
viruses (HRVs) that confer resistance to treatment with rupintri- 
vir, including the T129/130, T131, T143, N165, N130, and L136 
mutations (39). By checking the crystal structure of EV71 3CP™° 
(40), the N69S substitution unexpectedly appears to be close to 
but outside the active site of EV71 3CP™° (Fig. 8), which differs 
from the sites of the rupintrivir resistance-conferring mutations 
in HRV. The side chain of the N69 residue contacts with nearby 
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residues E71 and T132 via hydrogen bonds and stabilizes the con- 
formation of the active site of EV71 3C?*° when the substrate 
peptide or AG7088 inhibitor binds to 3CP"° (Fig. 8B and C). The 
replacement of the residue at this position by a serine residue 
shortens the side chain and would abolish these contacts (Fig. 8D). 
Interestingly, previous crystallographic work has found that the 
side chain of N69 forms an essential hydrogen bond with catalytic 
residue E71 to stabilize the catalytic triad, and biochemical analy- 
sis demonstrated that the interruption of this hydrogen bond by 
an N-to-S replacement significantly reduces the K,,,/k.,, values of 


m 


EV71 3CP™° (40). We therefore propose that this amino acid sub- 
stitution accommodates the higher structural flexibility of the ac- 
tive site of 3CP"° and allows polypeptide binding in the 3C”® active 
site in the presence of the inhibitor NK-1.8k. Since the residues 
constituting the active site of 3CP™° and the N69 residue are strictly 
conserved between EV71 and EV68 strains (Fig. 8E), it is conceiv- 
able that NK-1.8k has a consistent mechanism of inhibition of 
these two viruses. Determination of the complex crystal structures 
of NK-1.8k and WT 3CP™ and of the natural substrate and N69S 
3CP"° is warranted to dissect the precise role of N69 in the function 
of 30P™°, 
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